Decay of Loschmidt Echo Enhanced by Quantum Criticality 
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We study the transition of a quantum system S from a pure state to a mixed one, which is induced 
by the quantum criticality of the surrounding system E coupled to it. To characterize this transition 
quantitatively, we carefully examine the behavior of the Loschmidt echo (LE) of E modelled as an 
Ising model in a transverse field, which behaves as a measuring apparatus in quantum measurement. 
It is found that the quantum critical behavior of E strongly affects its capability of enhancing the 
decay of LE: near the critical value of the transverse field entailing the happening of quantum phase 
transition, the off-diagonal elements of the reduced density matrix describing S vanish sharply. 
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Introduction: Nowadays quantum -classical transitions 
described by a reduction from a pure state to a mixture 
P, 0] renew interests in many areas of physics, mainly 
due to the importance of quantum measurement and dc- 
coherence problem in quantum computing. To study this 
transition, some exactly-solvable models wereproposed 
for a system coupled to the macroscopic 0, 0, El or clas- 
sical || surrounding systems. Relevantly, in associa- 
tion with the quantum -classical transition in quantum 
chaos, the concept of Loschmidt echo (LE) from NMR 
experiments was introduced to describe the hypersensi- 
tivity of the time evolution to the perturbations expe- 
rienced by the surrounding system^, In this letter, 
by a concrete example, we will show how quantum phase 
transition (QPT) [lfj of the surrounding system can also 
sensitively affect the decay of its own LE, which means a 
dynamic reduction of its coupled system from pure state 
to a mixed one. Here, we note that a QPT effect has 
been explored for the Dicke model at the transition from 
quasi-integrable to quantum chaotic phases 

As a quantum critical phenomenon, QPT happens at 
zero temperature, at which the thermal fluctuations van- 
ish. Thus QPT is driven only by quantum fluctuation, 
and the uncertainty relation lie at the heart of various 
QPT phenomena. On the other hand, the randomness of 
the relative phase, which causes pure-mixed state tran- 
sition, also has its source in the uncertainty principle 
[l2j |. It is this observation that enlightens us to explore 
the relationship between QPT and the pure-mixed state 
transition described dynamically by the time evolution of 
the LE. It is common that the ground state of the critical 
system is very sensitive to the varying magnitude of the 
coupling constant or the system experiences a spon- 
taneous symmetry breaking at the critical point. Up to 
now, all of the known models of QPT possess this prop- 
erty. Actually, this kind of critical sensitiveness can be 
well understood by resorting to the concepts of quantum 
chaos through the LE Q or macroscopic enhancement of 
phase randomness [l2j |. 



Our approach is based on the Hepp-Coleman (HC) 
model [Si L| , which was initially proposed as a model for 
quantum measurement. In our generalization, the free 
spin 1/2 ensemble, as a model of measuring apparatus, is 
reDlaced bv the Isine SDin chain E in a transverse field. 

□ □ 




FIG. 1: (Color online) A schematic diagram of the physical 
implementation of the generalized Hepp-Coleman Model. The 
spins are arranged in a circle to form a ring array E. The 
central two-level system S possesses homogeneous couplings 
due to the overlaps of symmetric spacial wave function of S 
with those of spins. 

Model setup based on quantum phase transition: Our 
model concerning the decay of LE, illustrated in Fig. 1, 
is very similar to the Hepp-Coleman model 0, 0] or its 
generalizations We take the surrounding system 

E to be an Ising spin chain in a transverse field, which 
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satisfies the Born- Von Karman condition automatically. 
Consider a two level system S with the excited state |e) 
and the ground state \g), which is transversely coupled 
to E. The corresponding Hamiltonian reads as follows: 

H(X, 6)=-jJ2 foM+i + A<7| + S \e) (e| of) , (1) 



where J and A characterize the strengths of the Ising 
interaction and the coupling to the transverse field re- 
spectively; 6 indicates the small perturbation coupling of 
S to E; of (a = x, y, z) are the Pauli operators defined 
on the ith site in the lattice with lattice spacing a. 

We assume the two-level system is initially in a super- 
position state \4> s (0)) = c g \g) + c e |e), where the coeffi- 
cients c g and c e satisfy \c g \ 2 + |c e | 2 = 1. Then the evolu- 
tion of the Ising spin chain initially prepared in \ip (0)), 
will split into two branches \tp a (t)) — exp(—iH a t) \tp (0)) 
(a = e,g), and the total wave function is obtained as 
\ip (t)) = c g \g) <g> \tp g (t)) + c e |e) g> \ip e (t)). Here, the 
evolutions of the two branch wave functions \<p a (t)} are 
driven respectively by the two effective Hamiltonians 
H g = H(X,0) and H e = H(X,6) = H g + V e . Obviously, 
both H g and H e describe the Ising model in a transverse 
field, but with a tiny difference in the field strength. The 
central quantum system in two different states |e) and \g) 
will exert slightly different back actions on the surround- 
ing system, which manifest as two effective potentials 
V e = —J5J2j a j an d VslJt 0- This difference just results 
in the decay of the LE [9j defined as 

(2) 



L(X,t) = \{<p a {t)\<p e (*)> 



To show the key role of LE in depicting quantum deco- 
herence of the central system, and manifest the difference 
between LE and decoherence, we define the purity 0] 
P = Tr s (p 2 s ) = Tr s {[Tr E p(t)} 2 } to describe the de- 
coherence. Here, p(t) — (t)) (ijj (t)\, and Tr a means 
tracing over the variable of a, (a — E,S). A straight- 
forward calculation reveals the relationship between LE 
and the purity as P = 1 — 2 |c e c s | 2 [1 — L(X,t)]. This 
equation indicates that the purity depends on the initial 
states of the central system and the surrounding system 
E, but the LE only depends on the initial state of E. For 
simplicity, we assume the surrounding system is initially 
prepared in the ground state. In the following discussion, 
we will focus on the decay problem of LE induced by the 
coupling with the central system. 

Exact solution for the Loschmidt echo: We now pre- 
pare to prove that, just at the critical point A = A c = 1, 
the decay of LE is enhanced, accompanied by the QPT 
in one of the two evolution branches. 

We first diagonalize the effective Hamiltonian as H e = 
J2k £ e ( A k A k — 1/2^ in terms of the normal mode oper- 
ators ESQ 



^ — ikal , 



w*o 



] ), (3) 



which satisfy the canonical fermion anti-commutation 
relations. Here, TV is the number of sites of the spin 
chain, and op' = (— of±iof)/2 is defined by the 
Pauli matrices a", a — x,y,z. The coefficients u k e — 
cos (0g /2) , «g = sin (Q\ /2) depend on the the angle 



9 k e = 6 k e {5) = arctan[ 



— sin (ka) 
cos (ka) — (A + 6)~ 



(4) 



The corresponding single quasi-excitation energy e e is 

e k e (S) = 2J^J 1 + (A + Sf - 2 (A + 5) cos (ka). (5) 

Note that, in writing down the known result @ in a com- 
pact form, we have combined the Jordan- Wigner map 
and the Fourier transformation to the momentum space 

MM 

The effective Hamiltonian H g can be diagonalized in a 



similar way: H g = J2k £ g {j^l^k — 1/2J . In this case the 

single quasi-excitation energy is e k = £ k (0) and the cor- 
responding fermionic quasi-excitation operators can 
be obtained by the following Bogliubov transformation 



B± k = cos (a k ) A ±k - i sin (a k ) (A Tk 



(6) 



Here, a k = [0*-0*]/2, and 9 k are defined by 9 k = 8 k (0). 

We suppose that the spin chain is initially in the 
ground state \<f(0)) = \G) g of the Ising spin chain in 
a transverse field depicted by H g , i.e., B k \G) — for 
any operator B k . Then from Eq. © the state \G) can 
be rewritten as a BCS-like state: 



\G) g = 1] [coa(a k )-ism(a k )A{Al k ] \G) 



(7) 



k>0 



where \G) is the ground state of H e . This explicit ex- 
pression of \G) g enables us to calculate straightforwardly 
the LE (|2J|, which assumes the following factorized form: 

L(X, t) = Yl F k = - sin 2 (2a fc )sin 2 (e k t)}. (8) 



k>0 



k>0 



Quantum- classical transition at critical point of QPT: 
Since each factor F k in Eq (jHJ has a norm less than unity, 
we may well expect L(X,t) to decrease to zero in the 
large N limit under some reasonable conditions. This 
kind of factorized structure was first discovered and sys- 
tematically studied by one of the authors in developing 
the quantum measurement theory in classical or macro- 
scopic limit and it has been applied to analyze the 
universality of decoherence influence from environment 
on quantum computing |17| . Now we study in detail 
the critical behavior of the surrounding system near the 
critical point A c = 1 and its relation to the sensitive evo- 
lution of the LE perturbed by the central system even for 
a finite N. This turns out to reveal a novel mechanism 
responsible for the enhanced decay of LE. 
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Let us first make a heuristic analysis of the features of 
the LE. For a cut-off frequency K c we define the partial 
product for the LE 



L c (X,t) = l[F k >L(X,t), 



(9) 



fe>0 



and the corresponding partial sum S(X,t) — lnL c 



Unfit 



For small k we have e k e ft 



2~2k>0 

sin 2 [2a fc ] » ((5fca) 2 /(l - A) 2 (1 - A - <5) 
if c is small enough we have 



2J|1 — A — <5 1 , 
As a result, if 



S(A,t) 



S 2 E(K C ) sin 2 (2Ji|l -A-5|) 



(10) 



where E(K C ) = 4tt 2 N c (N c + 1)(2N C + 1)/(6N 2 ) and N c 
is the integer nearest to NK c a/2TT. Here we have used 
the fact that the Bloch wave vector k takes the discrete 
values 2mr/Na (n = 1, 2, •■ - N/2). In this case, it then 
follows that for a fixed t. 



L c (A,i)«exp(- 7 i 2 ) 



(11) 



when A — > A c = 1, where 7 = U 2 S 2 E(K C ) / '(1 - A) 2 . 

Notice that L(X,t) is less than L c (X,t). So from the 
above heuristic analysis we may expect that, when N 
is large enough and A is adjusted to the vicinity of the 
critical point A c — 1, the LE will exceptionally vanish 
with time. On the other hand, we observe that 7 seems 
to approach zero in the thermodynamic limit N — > 00 
for Na keeps as a constant in the process of taking this 
limit and E(K C ) cx 1/N 2 . Since a true QPT can occur 
just in the thermodynamic limit, it is natural to doubt 
whether the QPT, and thus the induced decay of the 
LE, can happen at the critical point. In fact, due to the 
vanishing denominator (1 — A) 2 of 7 at the critical point 
of the QPT, the decay of the LE is still possible even for 
7 having a vanishing numerator. For a practical system 
used to demonstrate the QPT induced decay of the LE, 
the particle number N is large, but finite, and then the 
practical 7 does not vanish. 

Now we resort to numerical calculation to test the 
heuristic analysis. For N = 50 ~ 250, 5 — 0.1, the LE are 
calculated numerically from the exact expression (JSJ) with 
the parameters within the ranges A S [0, 2], t G [0, 27/ J}. 
The results are demonstrated in Figs. 2a and 2b. 

In Fig. 2a there exists a deep valley in the domain 
around the line A = A c — 6 = 0.9. This reflects the fact 
that near the critical point of the surrounding system 
the LE is very sensitive to the perturbation experienced 
by the surrounding system. At the critical point, with a 
highly enhanced decay of LE, the central system transits 
from a pure state to a mixed state due to its entanglement 
with the surrounding system. The five curves in Figure 
2b clearly demonstrate the influence of N on the decay 
behavior of the LE. At A = A c — 6 = 0.9, the LE decays 
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FIG. 2: (Color online) (a) Three dimensional (3-D) diagram 
of the LE \L (A, t) 2 as the function of A and t for the system 
with N — 200. The valley around the critical point A c =l 
indicates that the decay of LE is enhanced by the QPT . The 
profile at A = 1 is in agreement with the analytical analysis, 
(b) The cross sections of the 3-D surface for the systems of 
N = 50, 100, 150, 200, and 250 at A= A c -5=0.9. It shows that 
the quasi-period of the LE is proportional to the size of the 
surrounding system . 



and revives as time increases. The period of the revival 
of the LE is proportional to the size of the surrounding 
system. 

Decays and revivals of Loschmidt echo as a witness of 
QPT: The novel phenomenon of the synchronization be- 
tween the QPT and the enhanced decay of the LE men- 
tioned above and its physical implication deserves fur- 
ther exploring. Generally, the two terms in H e represent 
two competitive physical effects with different order ten- 
dencies: in the weak coupling case A -C 1 the ground 
state is either all spins up or all spins down, while in the 
strong coupling case A 3> 1 the ground state tends to the 
saturated ferromagnetic state with all the spins pointing 
right. When A takes the value of the order unity, the 
qualitative properties of the ground states for A > 1 and 
A < 1 are similar to those for A ^ 1 and A <C 1 re- 
spectively. Only the critical point A — 1 has genuinely 
different properties. 

The singular behavior of QPT at A = A c reflects the hy- 
persensitivity of the ground states of the surrounding sys- 
tem with respect to the perturbation coupling imposed 
by the central system, which is reflected by the evolu- 
tion of LE. We can thus expect quantum evolution of the 
surrounding system to inherit this sensitivity, which can 
also be understood as a signature of quantum chaos: For 
a quantum system prepared in the identical initial state, 
two slightly different interactions can lead to two quite 
different quantum evolutions. Mathematically speaking, 
this means the LE, initially equal to 1, will decay with 
time and finally vanish. In this sense the sensitivity of 
quantum evolution to perturbation plays a crucial role 
in inducing the decay of LE. Due to the perturbations of 
two effective potentials caused by \e) and \g) respectively, 
the LE can decrease to zero due to the singularity at the 
critical point and the macroscopic enhancement of phase 
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FIG. 3: (Color online) The quantum phase transition en- 
hanced decay of LE at large N limit for small 5: except for 
5 = 0.01 the explanations are the same as that in Fig. 2. 

randomness for large N, only at which QPT occurs [l^ . 

Next, we numerically investigate the LE for system 
with finite N. It turns out that as N increases the LE 
will vanish for very small 6. For example, we take S = 
0.01, N = 500 ~ 2500 and compare the numerical results 
illustrated in Fig. 3 with those for S = 0.1, N = 50 ~ 250 
in Fig. 2. From Fig. 2a and 3a one can clearly see that 
the valley narrows as 5 decreases and N increases. On 
the other hand, the comparison between Fig. 2b and 
3b shows an interesting phenomenon in the LE at the 
critical point A c . Firstly, the periods of the revival of 
the LE for the two cases are both proportional to the 
size of E. Secondly, the behaviors of the LE in Fig. 2b 
and 3b are almost same. The numerical results seem 
to suggest an extrapolation for a scaling behavior, i.e., 
the LE at the critical point L c (t, 5, N) is invariant under 
the scaling transformation t — > t/a, 8 — > aS, and N — > 
N/a. This prediction is in agreement with the analytical 
analysis with the small Jt approximations. Therefore, 
our result opens a possibility that the distinctive picture 
of the decays and revivals of the LE may serve as a good 
witness of QPT in the case of finite N. 

Conclusion: In summary, by a special model, we have 
obtained the exact expression of the LE and analyzed 
the possible relation between the quantum-classical tran- 
sition of the central system, characterized by the LE of 
the surrounding system, and the occurrence of a QPT 
in its surrounding system E. Both the heuristic analy- 
sis and the numerical calculations we performed reveal a 
novel mechanism of the decay of the LE. It is well known 
that quantum critical phenomenon is closely associated 
with the entanglement among the qubits consisting of the 
surrounding system . We would like to stress that our 
present study is from a different perspective; our empha- 
sis is on the relation between the QPT of the surrounding 
system and its entanglement with the central system S, 
which is qualitatively characterized by the LE. In our 
model, the maximal quantum entanglement between S 
and E can be reached and then the central system tran- 
sits from a pure state to a mixed state when a QPT 



of E takes place in one of the two evolution branches. 
This result suggests an unexplored and rather intriguing 
relationship among entanglement, LE, decoherence and 
criticality. 
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